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Abstract. Population protocols have been introduced as a model of 
sensor networks consisting of very limited mobile agents with no control 
over their own movement. A population protocol corresponds to a col- 
lection of anonymous agents, modeled by finite automata, that interact 
with one another to carry out computations, by updating their states, 
using some rules. 

Their computational power has been investigated under several hypothe- 
ses but always when restricted to finite size populations. In particular, 
predicates stably computable in the original model have been character- 
ized as those definable in Presburger arithmetic. 

We study mathematically the convergence of population protocols when 
the size of the population goes to infinity. We do so by giving general 
results, that we illustrate through the example of a particular population 
protocol for which we even obtain an asymptotic development. 
This example shows in particular that these protocols seem to have a 
rather different computational power when a huge population hypothesis 
is considered. 



1 Motivation 



The computational power of networks of finitely many anonymous resource- 
limited mobile agents has been investigated in several recent papers. In partic- 
ular, the population protocol model, introduced in pQ, consists of a population 
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of finite-state agents that interact in pairs, where each interaction updates the 
state of both participants according to a transition based on the previous states 
of the participants. When all agents converge after some finite time to a common 
value, this value represents the result of the computation. 

According to survey [3J, the model corresponds to the following assump- 
tions: anonymous finite-state agents (the system consists of a large population 
of indistinguishable finite-state agents), computations by direct interactions (an 
interaction between two agents updates their states according to a joint tran- 
sition table), unpredictable interaction patterns (the choice of interactions is 
made by an adversary, possibly limited to pairing only agents in an interaction 
graph), distributed input and outputs (the input to a population protocol is 
distributed across the initial state of the entire population, similarly the output 
is distributed to all agents), convergence rather than termination (the agent's 
output are required to converge after some time to a common correct value). 

Typically, in the spirit of [1] and following papers (see again [3] for a sur- 
vey), population protocols are assumed to (stably) compute predicates: a pop- 
ulation protocol stably computes a predicate <f>, if for any possible input x of 
<j>, whenever 4>{x) is true all agents of the population eventually stabilize to a 
state corresponding to f , and whenever <fi(x) is false, all agents of the population 
eventually stabilize to a state corresponding to 0. 

Predicates stably computable by population protocols in this sense have been 
characterized as being precisely the semi-linear predicates, that is to say those 
predicates on counts of input agents definable in first-order Presburger arithmetic 
[TTj . Semilinearity was shown to be sufficient in [1] and necessary in [2] . 

Here, we study a new variant: we assume the number of agents in the popu- 
lation to be huge, or close to infinity (we call this a huge population hypothesis), 
and we don't want to focus on protocols as predicate recognizers, but as comput- 
ing functions. We assume outputs to correspond to proportions, which are clearly 
the analog of counts whenever the population is infinite or close to infinity. 

We do so by providing general results that we illustrate by considering a 
particular population protocol, that we prove to converge to a fraction of 
agents in a given state whatever its initial state is. We hence show that some 
algebraic irrational values can be computed in this sense. Then we show how 
the reasoning behind the proof of convergence of this particular protocol can be 
generalized to any such protocol to prove that its behaviour can be analyzed 
through use of deterministic differential equations. We also give an asymptotic 
development of the convergence in the case of the protocol computing 

Our motivation is twofold. First, to prove formally that population protocol 
with a huge population hypothesis can be mathematically studied using popula- 
tion models and ordinary differential equations. Second, to show that protocols 
considered with these two hypotheses (huge population, computing functions 
and not only predicates), have a rather different power. 

We consider this work as a first step towards understanding which numbers 
can be computed by such protocols. Whereas we prove that can be computed, 
and whereas this is easy to see that computable numbers in this sense must 



be algebraic numbers of [0,1], we didn't succeed yet to characterize precisely 
computable numbers. 

In this more long term objective, the aim of this current work is first to 
discuss in which sense one can say that these protocols compute an irrational 
algebraic value such as and second to study mathematically formally the 
convergence. 

Our discussion is organized as follows. In Section [51 we discuss related work. 
In Section[3l we recall the model of population protocols. In Section[4j we present 
a particular population protocol and we explain in which sense we would like 
to say that this protocol computes some irrational algebraic value, with a huge 
population hypothesis. We do so first by some informal study, that we justify 
mathematically in the following three sections. We use a general theorem about 
approximation of diffusions presented in Section [3] How this theorem yields the 
proof of convergence is explained in Section [6] Then, in Section [3 we show how 
the method used on this particular example can be extended to prove that the 
study of such protocols reduces to the study of differential systems. The following 
two sections deal with additional results that can be given on restricted parts 
of the model. We prove in Section [H] that this is even possible to use the same 
theorem to go further and get an asymptotic development of the convergence 
on the example of Section [9] gives, without proof, some first results on the 
type of algebraic numbers that can be computed in this model using only two 
possible states for all agents. Finally, Section [TO] is devoted to a conclusion and 
a discussion. 

2 Related Work 

The population protocol model was inspired in part by the work by Diamadi 
and Fischer on trust propagation in social networks [7] . The model is introduced 
and motivated in [T] by the study of sensor networks in which passive agents are 
carried along by other entities. The canonical example given in this latter paper 
is the example of sensors attached to a flock of birds. 

So far, most of the works on population protocols have concentrated on char- 
acterizing what predicates on the input configurations can be stably computed in 
different variants of the models and under various assumptions. In particular, in 
PQ, Angluin et al. proved that all semi-linear predicates can be computed by the 
original model but left open the question of their exact power. This was solved 
in [3J, where it has been proved that no- more predicates can be computed. 

Variants of the original model considered so far include restriction to one- 
way communications, restriction to particular interaction graphs, random in- 
teractions, self-stabilizing solutions, through population protocols, to classical 
problems in distributed algorithmic, the analysis of various kind of failures of 
agents an d their consequences, etc. See survey [3]. 

As far as we know, a huge population hypothesis in the sense of this pa- 
per, has not been considered yet, except in very recent papers |4|5j . In these 



papers, the authors also study the dynamics and stability of probabilistic popu- 
lation protocols via an ordinary differential equation approach. They introduce 
several continuous models inspired by population models in mathematics and 
by population protocols. Some relations between continuous models are proved. 
Results from theory of stability for continuous systems are used to discuss the 
stability of protocols. The models considered in [415] are not provably connected, 
in these papers, to classical finite population protocols. We prove in this paper 
that there is a strong relation between classical finite population protocols and 
models given by ordinary differential equations, and hence with the models of 
4 5 . We also prove that the macroscopic view given by ordinary differential 
equations and by population models is provably the weak limit of classical finite 
population protocols when the size of the population goes to infinity. Compared 
with |4l5j . we focus on computability by these models, whereas |4I5] focus on 
stability questions, and do not discuss computability issues. 

Notice that we assume that interactions happen in a probabilistic way, ac- 
cording to some uniform law. In the original population protocol model, only 
specific fairness hypotheses were assumed on possible adversaries pQ. When the 
size of the population goes to infinity, uniform sampling of agents seems to us the 
most natural way to extend the fairness hypothesis. This assumption is consis- 
tent with the interpretation of agents as autonomous biological entities moving 
at random. Notice that this notion of adversary has already been investigated 
for finite state systems [3]. 

The results proved in this paper can be considered as a macroscopic ab- 
straction of a system given by microscopic rules of evolutions. See survey [3] for 
general discussions about extraction of macroscopic dynamics. 

Whereas the ordinary differential equation ((5]) can be immediately abstracted 
in a physicist approach from the dynamic ([!]) , the formal mathematical equiva- 
lence of the two approaches is not so trivial, and is somehow a strong motivation 
of this work. 

Actually, these problems seem to arise in many macroscopic justifications 
of models from their microscopic description in experimental science. See for 
example the very instructive discussion in [10] about assumptions required for 
the justification of the Lotka-Volterra (predator-prey) model of population dy- 
namics. In particular, observe that microscopic correlations must be neglected 
(i.e. £[J7] = is needed, where E denotes the expectation). With a 

rather similar hypothesis (here assuming E[p 2 ] = E[p] 2 ), dynamic ([5]) is clear 
from rules ([1]). Somehow, we prove here that this last hypothesis is not necessary 
for our system. 

The techniques used in this paper are based on weak convergence techniques, 
introduced in [12] , relating a stochastic differential equation (whose solutions 
are called diffusions) to approximations by a family of Markov processes. More 
specifically, we use a weaker form [6] of a theorem from [12] . Refer also to [8] for 
an introduction to these techniques. 



3 Population Protocols 



We now recall definitions from pQ. A protocol is given by (Q, S, l, u, S) with 
the following components. Q is a finite set of states. S is a finite set of input 
symbols, i : S — > Q is the initial state mapping, and uj : Q — ► {0, 1} is the 
individual output function. 6 C Q 4 is a joint transition relation that describes 
how pairs of agents can interact. Relation <5 is sometimes described by listing all 
possible interactions using the notation ((71,92) — > {^11^2)1 or even the notation 
9i<72 — > 9W2' f° r (9i j 92, ei, 9 2 ) e ^ (with the convention that (91,92) — > (91,92) 
when no rule is specified with (91, 92) in the left hand side). 

Computations of a protocol proceed in the following way. The computation 
takes place among n agents, where n > 2. A configuration of the system can 
be described by a vector of all the agent's states. The state of each agent is an 
element of Q. Because agents with the same states are indistinguishable, each 
configuration can be summarized as an unordered multiset of states, and hence 
of elements of Q. 

Each agent is given initially some input value from S: Each agent's initial 
state is determined by applying 1 to its input value. This determines the initial 
configuration of the population. 

An execution of a protocol proceeds from the initial configuration by inter- 
actions between pairs of agents. Suppose that two agents in state 91 and 92 meet 
and have an interaction. They can change into state q[ and q' 2 if (91, 92, q[, q'2) is 
in the transition relation 8. If C and C are two configurations, we write C — > C 
if C can be obtained from C by a single interaction of two agents: this means 
that C contains two states 91 and 92 and C is obtained by replacing 91 and 92 
by q[ and q' 2 in C, where (91,92,91,92) € 5. An execution of the protocol is an 
infinite sequence of configurations Co, C±, C%, ■ ■ ■ , where Co is an initial configu- 
ration and Ci — > Cj-i-i for all i > 0. An execution is fair if for all configurations C 
that appears infinitely often in the execution, if C — > C for some configuration 
C, then C appears infinitely often in the execution. 

At any point during an execution, each agent's state determines its output 
at that time. If the agent is in state 9, its output value is oj(q). The configuration 
output is (respectively 1) if all the individual outputs are (respectively 1). If 
the individual outputs are mixed 0s and Is then the output of the configuration 
is undefined. 

Let p be a predicate over multisets of elements of S. Predicate p can be con- 
sidered as a function whose range is {0, 1} and whose domain is the collection 
of these multisets. The predicate is said to be computed by the protocol if, for 
every multiset /, and every fair execution that starts from the initial configura- 
tion corresponding to J, the output value of every agent eventually stabilizes to 
P(I)- 

The following was proved in |H2j 

Theorem 1 ( [T1I2] ) . A predicate is computable in the population protocol model 
if and only if it is semilinear. 



Recall that semilinear sets are known to correspond to predicates on counts 
of input agents definable in first-order Presburger arithmetic |f lj . 



4 A Simple Example 

Consider the following population protocol, with Q = {+, — }, and the following 
joint transition relation. 



(1) 



Using previous (classical) definition, this protocol does not stably compute 
anything. Indeed, if we put aside the special configuration where all agents are 
in state — which is immediately left in any next round, any configuration is 
reachable from any configuration. 

However, suppose that we want to discuss the limit of the proportion p(k) 
of agents in state + in the population at discrete time k. If n+(fc) denotes the 
number of agents in state +, and n_(fc) = n — n + (k) the number of agents in 
state — , 

P (k) = !^*>. 

n 

From now on, we suppose that at each time step, two different agents are 
sampled uniformly among the n particles, independently from the past. Since we 
are dealing with n indistinguishable agents, the population protocol is completely 
described by the number of agents in state +. We are then reduced to determine 
the evolution of the Markov chain 

The above discussion ensures that (p(k)) is an irreducible Markov chain in 
{—,...,—}. Let us now compute the transition probabilities of this irreducible 
Markov chain. We have 

p(fc + i)-p(fc)e{-i,i}. 

Then, we have to determine for each i = 1,2, ... ,n 

* (f °(£ - ^) ■■= P (P(k + 1) = ^ I P(k) = % - 

* (n) -> l -ir) ■= p (p( k + 1) = — I P( k ) = ~ 
\ n n 

Assume that p(k) = i/n : p(k) decreases only if the two agents sampled are in 
state +. That is, 



JT (n)(± = ill = ~ ^ 



G) n(n-l)' 



In any other case, p(k) increases by one : 

7i-(")f A — > = l - 7r( n )f-i. _> i=l) 

n(n — 1) ' 

A consequence of the ergodic theorem is that the chain (p(fc)) admits a unique 
stationary distribution /U. By definition, it is the only application 

-}-[0,l] 

such that 
!• E"=i /*(*/") = L 

2. /x satisfies the balance equation, i.e. for each i 

n n n n n 

We do not pay attention to the exact expression of fi. We only notice that, as 
the unique solution to a rational system, it is an clement of Q™. Hence, its mean 
n{i/n)i/n is a rational number, that we denote p^ n \ 

A second consequence of the ergodic theorem is the following convergence : 

pW+pW + -+pW fe ^V«>, almost surely. 
k 

The purpose of the rest of the discussion is to show that however, when n goes 
to infinity, the mean value of p{k) converges to the irrational number V2/2. To 
see why this can be expected, observe that 

V " n ' F K ' n-1 iK n- 1' 

and write 

E[n + (k + 1) - n+(k) | n+(fc)] = 7rW(i - ^) - i=l) 

= l_ 27r («)(A^i^) (2) 

= l-V(fc)^r+Mfc)^T 

From this, we can derive (not yet rigourously) the asymptotic behavior of 
p(k). Take indeed n large, so that the right-hand term is close to 1 — 2p(k) 2 . Now, 
when k goes large, if there is some convergence of the mean proportion of +, the 
system must concentrate on configurations that does not create or destroy +, in 
mean. Thus, 1 — 2p(k) 2 must vanish and p(k) w V2/2. 

Hence the remaining problem is to justify and discuss mathematically the 
convergence. 



5 A General Theorem about Approximation of Diffusions 

We will use the following theorem from . We use here the formulation of it 
in [Sj (Theorem 5.8 page 96). 

Suppose that for all integers n > 1, we have an homogeneous Markov chain 
(Xk ) m ^ d w ith transition kernel 7r'™' (x, dy), meaning that the law of Y^f^, 
conditioned on Yq , • • • , K. , depends only on K. and is given, for all Borelian 
B, by 

P(Y^\eB\Y^) = ^(Y^,B), 

almost surely. 

Define for x € K d , 



6 (n) (z) 


= "/ 


f (y-x)n^(x 1 dy), 


a (n) (x) 


-J 


(y~ x){y- x)*n {n) (x,dy), 


^(z) 


-J 


r (y-x) 3 7r^(x,dy), 




= nir { 


' n \x,B{x,e)% 



where B(x, e)° denotes the complement of the ball with radius e, centered at x. 
In other words, 

b^(x) = nE x [(Y 1 -x)], 

and 

a^\x)=nE x [{Y l -x){Y l -xf] 
where E x stands for "expectation starting from x" , that is, 

E x [{Y 1 -x)]=E[(Y 1 -x)\Y = x}. 

The coefficients and can be interpreted as the instantaneous drift and 
the variance (or matrix of covariance) of X^ n \ 
Define 

* {n) (t) = y$ 1 + N- lnt])(Y^ +11 -yW). 

Theorem 2 (Theorem 5.8, page 96 of |6j). Suppose that there exist some 
continuous functions a, b, such that for all R < +oo, 

lim supi x |< fl |a (n) (a;) - a(x)\ = 

n — ►do 

lim sup lxl<R \b {n) (x) - b(x)\ = 
n— >oo 

lim «up w<H 2lW = 0,Ve > 

n — >oo 

sup K^ n \x) < oo. 

|x|<iJ 



With a a matrix such that o~(x)a* (x) — a(x), x € M. d , we suppose that the 
stochastic differential equation 

dX{t)=b{X{t))dt + (j{X{t))dB{t), X(0)=x, (3) 

has a unique weak solution for all x. This is in particular the case, if it admits 
a unique strong solution. 

Then for all sequences of initial conditions Y^ — > x, the sequence of random 
processes converges in law to the diffusion given by . In other words, for 
all function F : C(R + ,R) — > K bounded and continuous, one has 

lim E[F(X^)} =E[F(X)}. 

n — >oo 

6 Proving Convergence of Previous Simple Example 

Consider Y^ as the homogeneous Markov chain corresponding to p(k), when n 
is fixed. From previous discussions, (x, .) is a weighted sum of two Dirac that 
weight x — — and x + — , with respective probabilities tt_i and 7r +1 , whenever a; 
is of type j % for some i. 

From Equation @ we have 

1 T7 9 

+ 1) -p(k)\p(k)} = -(1 - 2p{kf -+p(k) -), (4) 

n n — 1 n — 1 

which yields 

= 1 - 1p(kf-^- + P (k)^—, 
n — 1 ri — l 

when x = i/n. Now, clearly (p(k + 1) — p(k)) 2 = \, and hence 

a< B >(aO = -, 

71 

when x = i/n. Taking a(x) = and b(x) = 1 — 2x 2 , we get 
lim sup\ x \<n\a( n ' (x) — a(x)\ = 

lim sup\ x \ <R \b^ n \x) - b(x)\ = 

for all R < +00. Since the jumps of are bounded in absolute value by — , 
is null, as soon as - is smaller than e, and so 

lim sup lxl<R A( n) = 0,Ve > 0. 

Tl-*CO 

Finally, 

sup < oo 

\x\<R 



is easy to establish. 

Now, (ordinary and deterministic) differential equation 

dX(t) = (1 - 2X 2 )dt (5) 

has a unique solution for any initial condition. 

ft follows from above theorem that the sequence of random processes X( n > 
defined by 

* (n) w=n3 +(»*- Ln*j)(yjS2fij - Y Si) 

converges in law to the unique solution of differential equation (|SJ • 

Clearly, all solutions of ordinary differential equation $5$ converge to 
Doing the change of variable Z(t) — X(t) — we get 

dZ(t) = (-2Z 2 + 2V2Z)dt, (6) 

that converges to 0. 

Coming back to p(k) using definition of X {n \t), we hence get 

Theorem 3. We have for all t, 

P([nt\) = \ + Z n (t), 

where Z n (t) converges in law when n goes to infinity to the (deterministic) solu- 
tion of ordinary differential §§§ . Solutions of this ordinary differential equation 
go to at infinity. 

This implies that p(k) must converge to ^ when k and n go to infinity. 



7 Generalization To General Population Protocols 

We will now generalize the reasoning made on this particular example in order 
to prove that the behaviour of any population protocol can be approximated by 
a deterministic differential equation in a similar way. 

Transition rules of a population protocol are of the form: 

q q' ^ S 1 (q,q / ) S 2 (q,q') 

for all (qi,g 2 ) G Q 2 . 

As previously we consider pairwise interactions between two agents chosen 
randomly according to an uniform law in a population of size n. 

Let us define the Markov chain corresponding to the vector of whose 
components are the proportions of agents in the different states and 

^ n \t)=Y^+(nt-[nt\)(Yf:l 1 ,-Y{:h. 



Theorem 4. Le b be the function defined by : 

b ( X )= X q X q'(-( e q+ e 'q)+ e Sl(q.,q')+ e S 2 (q..q')) 

{q,q r KQ 

where (e g ) g6 Q is the canonical base of MP. 

Then for all sequences of initial conditions Yq — > X, the sequence of ran- 
dom processes XW converges in law to the solution of the ordinary differential 
equation: 

dX(t) = b(X(t))dt, X(Q) = x, (7) 

Remark 1. Ordinary differential equation |(7J| corresponds to a degenerated stochas- 
tic differential equation. Being deterministic, we are sure that it has a unique 
weak solution for all x. 

Proof. Y™ is of the form required by Theorem [5] with n^ n '(x, .) being the sum of 
5IQI Dirac : the variation of the proportion of agents in any given state belongs 
to { — , — , 0, —, — , } and the probabilities of any of these variations are clearly 
only dependant on the current state x. 

Now let us define a^ n \x),b^{x),K i - n \x) and A ( e n) as in Theorem H Let R 
be any finite non-negative real number. 

As in the example above, since at any given time step at most two out of n 
agents change state, A™ — if e > — and thus 

lim sup ]x \< r AW = 0,Ve > 0. 



sup K {n) {x) < 00 

\x\<R 

is also easy to establish. 
Similarly 

Vz G K IQI , \x\ < R, \a {n) (x)\ < 

n 

So if we take a(x) = 0, we have 

lim supu|<fi|a (n) (a;) - a(x)\ = 0. 

n — >oo 

If we write, for all (q, q') G Q 2 , q 7^ q', 

HqJ' ( X ) = x q x q' _ T 

and 



^■q.q( X ) ~ x q x 



n X n 



n — 1 n—1 



Then IJ^ n q , (x) is exactly the probability of an encounter between an agent in 
state q and an agent in state q' to happen when the population is in configuration 
x. We then have : 

&W(x)= ^W(-( e 9 + <) + ^(^')+^(^'))' 

(<?,<?')6<3 

or 

b (n \x) = ^ZjHx) - X l(~ 2e l + e Sl(q,q) + e <M<?,9))- 

qeQ 

Thus, finally, 

lim sup\ x \ <R \b [n) (x) - b{x)\ = 0. 

n — >oo 

We can now conclude by Theorem [2J 

This means that to understand the asymptotic behaviour of any such pro- 
tocol, we can study the associated differential equation, ft is also of interest to 
note that the function b defined here is a quadratic form over M9 . 



8 An Asymptotic Development of the Example Dynamic 

It is actually possible to go further, at least in some cases like the previous simple 
example, and prove the equivalent of a central limit theorem, or if one prefers, 
to do an asymptotic development of the convergence, in terms of stochastic 
processes. We shall do so for the example dynamic used before using the same 
notations as in previous sections. 

In our previous example, as p(k) is expected to converge to consider the 
following change of variable: 

Y^(k) = V^( P (k)-^-). 

The subtraction of is here to get something centered, and the y/n factor 
is here in analogy with classical central limit theorem. 

Clearly, Y^ n '(.), that we will also note Y(.) in what follows when n is fixed, 
is still an homogeneous Markov Chain. 

We have 

E[Y(k + 1) - Y{k)\Y{k)] = V^(E[ P (k + 1) - p(k)\p(k)}), 
hence, from Q, 

1 n 2 

E[Y(k + 1) - Y(k)\Y(k)} = -=(1 - 2p{kf +p(k) -). 

Jn n — 1 n — 1 



Using = ^ + , we get 

E [Y(k + 1) - = + + m^tt) + ^w 2 (-^i) ) 

which yields the equivalent 

nE[Y(k + 1) - Y{k)\Y{k)\ w -2^2 F(fc) 
when n goes to infinity. We have 

£[(Y(fc + 1) - r(fc)) 2 |F(fc)] = n(E[(p(k + l)-p(k)) 2 \p(k)}), 
hence, since (p(k + 1) — p{k)) 2 is always 

n£[(F(fc + 1) - Y(k)) 2 \Y(k)} = 1. 
Set a(x) = -2y/2x, b(x) = 1. Fr om the above calculations we have clearly 
lim sup\ x \<:^\a^ n \x) — a(x)\ = 

n — >oo 

lim siipi.Kfll^'ta;) - b(x)\ = 

for all < +00. Since the jumps of are bounded in absolute value by -4=, 
is null, as soon as is smaller than e, and so 

lim «i P | x |< fl 4 n) =0,Ve>0 

n — >oo 

sup i^ (n) (x) < oo 

\x\<R 

is still easy to establish. 

Now the stochastic differential equation 

dX{t) = -2V2X(t)dt + dB(t) (8) 

is of a well-known type. Actually, for b > and u / 0, stochastic differential 
equations of type 

dX{t) = -bX(t)dt + adB(t) 

are known to have a unique solution for all initial conditions X(0) = x. This 
solution is an Orstein-Uhlenbeck process, and is given by (see e.g. [6]) 

X(t) = e- bt X{0) + [ e-^-^odBis). 
Jo 

For all initial conditions X(0), X(t) is known to converge in law, when t 
goes to infinity, to the Gaussian distribution Af(0, |r). This latter Gaussian 
distribution is an invariant distribution for the Orstein-Uhlenbeck process. See 
for example [6]. 

We have all the ingredients to apply Theorem [2] again. We obtain: 



Theorem 5. For all t, 

[2 1 
p{\jit\) = \ + - r A n {t), 
2 y/n 

where (A„ (t)) t >o converges in law to the unique solution of stochastic differential 
equation ([8]). Hence A n (t) converges toM(0, ^) when t and n go to infinity. 



9 Some Other Algebraic Numbers 

We have treated in detail the case of y/2/2. We present in this Section, without 
proofs, the extension of our result to 2-states protocols, with pairing of agents. 

At each time step, two agents are picked and their states are possibly changed, 
according to the fixed rule S. These protocols are completely described by the 
three mean increments 

a = E[n + (k + 1) — n + (k)\{+, +}have been picked], 
(3 = E[n + (k + 1) — n+(fc)|{+, — }have been picked], 
7 = E[n + (k + 1) — n + (k)\{— , — }have been picked]. 

Thus, there are 27 = 3 3 different rules. We denote them by the corresponding 
triplet (a, [3, 7). For instance, the rule computing V2/2 is denoted by (— f , +1, +1). 
We exclude the identity rule (0,0,0). 
We also set 

a = a — 2/3 + 7, 
b = 2(3- 2 7 , 
c = 7. 

Then, we associate to each triplet (a, (3, 7) the polynomial 

P = aX 2 + bX + c. 

The following Lemma, whose proof is omitted, is the basis of the next discussion. 



Lemma 1 The polynomial P = aX 2 + bX + c admits at most one root in (0, 1), 
which we denote by p* . Moreover, 

q := (a 2 - 2(3 2 + 7 2 )(;P*) 2 + (2/? 2 - 2 7 V + 7 2 > 

and 

2ap* + b<0. 



We will see that the computational power of a protocol population reads on the 
corresponding polynomial P. 



Case 1: P has no root in (0,1). Monotonic convergence. 

For 10 rules, P does not admit a root in (0, 1). In this case, the convergence 
of the corresponding protocol is easy to establish. Take for instance (0, 1,2): 

'++ + 
< +- ^ + + 
~+ + 

It is clear that the protocol converges to the configuration {+}". We summarize 
the behaviors of the 9 remaining rules in the following table. 



a 


a 


7 


Convergence 





i 





{+}" ( or {~ }" ^ ^ i s the init. config.) 





i 


1 


{+}" 





i 


2 


{+}" 








1 


{-H+}" -1 








2 


{+r or {-H+}"- 1 





-1 





{— } n (or {+}" if it is the init. config.) 

{+H-}"- 1 


-1 








-1 


-1 





{-}" 


-2 








{-r ou i+ji-r- 1 


-2 


-1 








Case 2: P has a unique root in (0,1). Approximation with a diffusion. 

Depending on (a, [3, 7), p* has one of the three following expressions: 

-b+Vb' 2 -4ac 

2a ' 

~* _ ) — b— Vfc 2 — 4ac 
V ] 2a ' 

c 

As in the previous section, we set 

n = n (n) ==^(pi n) -P*), 

and X is the linear interpolation of F: 



Theorem 6. Assume that p^ converges to a random variable X in (0,1). 
When n goes to infinity, the process (X^(t)) 4>Q converges to the unique (weak) 
solution X of the stochastic differential equation 

dX t = {2ap k + b)X t dt + qdB t (9) 

starting at Xq. Since, by the previous Lemma, 2ap* + b < and q > 0, this solu- 
tion X is an Ornstein- Uhlenbeck process, and has the following representation: 

X t = Xoe^+W + q f exp^+^-s) dB s . 
Jo 

Proof. We omit the proof, for it extends easily from the case y/2/2. 

There are 16 rules for which P has a root in (0, 1). They compute 13 different 
algebraic numbers, as shown in the next tabic. 



a 




7 


Polynomial P 


p* 





-1 


1 


3X 2 - AX 


+ 1 


1/3 





-1 


2 


AX 2 - 6X 


+ 2 


1/2 


-1 


1 





-3X 2 + 


2X 


2/3 


-1 


1 


1 


-2X 2 


+ 1 


a/2/2 


-1 


1 


2 


-X 2 - 2X 


+ 2 


V3 - 1 


-1 





1 


-2X 


+ 1 


1/2 


-1 





2 


X 2 -AX 


+ 2 


2-V2 


-1 


-1 


1 


2X 2 - AX 


+ 1 


1 - v/2/2 


-1 


-1 


2 


3X 2 - 6X 


+ 2 


1 - V3/3 


-2 


1 





-AX 2 + 


2X 


1/2 


-2 


1 


1 


-3X 2 


+ 1 


V3/3 


-2 


1 


2 


-2X 2 - 2X 


+ 2 


(V5-l)/2 


-2 





1 


-X 2 - 2X 


+ 1 


a/2-1 


-2 





2 


-AX 


+ 2 


1/2 


-2 


-1 


1 


X 2 -AX 


+ 1 


2- V3 


-2 


-1 


2 


2X 2 - 6X 


+ 2 


(3 - VE)/2 



10 Conclusion 

Population protocols have been introduced in 1J as a model for sensor net- 
works. In this paper we considered population protocols with a huge population 
hypothesis. 

Whereas for the definitions of computability considered in pQ, some popula- 
tion protocols are not considered as (stably) convergent, we proved through an 
example that sometimes they actually compute in some natural sense some irra- 
tional algebraic values: indeed, we gave a simple example where the proportion 
of agents in state + converges to whatever the initial state of the system is. 



One aim of this paper was to formalize the proof of convergence. We did it on 
this example and on general population protocols using a diffusion approximation 
technique, using a theorem due to [T2] . 

We consider this work as a first step towards understanding which numbers 
can be computed by such protocols in this sense. We gave some preliminary 
results about numbers computable with 2-states protocols. 

Considering more states, it is easy to derive from the protocol considered here 

another protocol that would compute *J by working with an alphabet made 

of pairs of states. Whereas it is easy to see that computable numbers in this 
sense must be algebraic numbers of [0, 1], wc didn't succeed yet to characterize 
precisely computable numbers. 
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